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1. Introduction 

The Virial Theorem (VT) has been known for a long time in both classical mechanics 
and quantum mechanics. In the classical case, it provides a general equation relating 
the average over time of the kinetic energy (T) with that of the function of potential 
energy (r ■ W). The VT was given its technical definition by Clausius in 1870 [1]. 
Mathematically, the theorem states 

2{T) = {f.VV). (1) 

If the potential takes the power function V{r) = ar^ with r = |r|, the VT adopts a 
simple form as 

2{T)=n{V). (2) 

Thus, twice the average kinetic energy equals n times the average potential energy. The 
VT in quantum mechanics has the same form as the classical one, except for the average 
over time in Eqs. ([T]) and ^ replacing by the average over an energy eigenstate of the 
system. It dates back to the old papers of Born, Heisenberg and Jordan [2], and is 
derived from the fact that the expectation value of the time- independent operator r ■ p 
under a eigenstate is a constant [3], 

^j^{^\r-m = {nr-p,H]\^) = 0, (3) 

2 

where H = ^ + V is the Hamiltonian and \ip) is an eigenket of H. 

In 1960, Hirschfelder [1] generalized the relationship by pointing out that f-p could 
be replaced by any other operators which were not dependent on time explicitly. In this 
way, he established the Hypervirial Theorem (HVT). For example, in a one-dimensional 
system, one can replace f-p = xp hj the hypervirial operator x^p, and obtain the 
recurrence relation of (x^), 

2kE{x''-^) = 2k{x^~W) + (x'^^) - hi{k - l){k - 2){x''-^), (4) 

where k is an integer and E is the eigenenergy. 

The Hellmann-Feynman (HF) Theorem is another important theorem in quantum 
mechanics, which has been applied to the force concept in molecules by using the 
internuclear distance as a parameter [5l E]. Let the Hamiltonian H{^) of a system 
be a time-independent operator that depends explicitly upon a continuous parameter 
^, and lipiO) be a normalized eigenfunction of H{^) with the eigenvalue i.e. 
HiOl^iO) = EmiOl^iO)^ (V'(0IV'(0) = 1- The HF theorem states that 

If the potential takes the power function V{r) = ar", the HF gives an equation 
representing the relation between eigenenergy E^ and mean value of r". 
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Figure 1. (a) The gnomonic projection, which is the projection onto the tangent 
plane from the center of the sphere in the embedding space, (b) Two classical orbits 
are shown as si (dashed line) and S2 (solid line), corresponding to the circular motion 
and the radial one respectively. 



Based on the relations in Eqs. ^ and ([6]), the Hypervirial-Hellmann-Feynman 
Theorem (HVHF) perturbation theorem is estabhshed [3 E]. It provides a very 
efficient algorithm for the generation of perturbation expansions to large order, replacing 
the formal manipulation of Fourier series expansions with recursion relations. This 
perturbation method just need the energy instead of the wave functions of the system, 
and it is easy to achieve on the computer. 

These results are well known, but have not, to our knowledge, been exploited in a 
curved space. In the present work, we focus on the one- and two-dimensional spherical 
geometry. The coordinate systems adopted in this paper are shown in Fig. [l](a): (i) An 
intuitive way to describe a two-dimensional sphere is to embed it in a three-dimensional 
Euclidean space. Each pair of independent variables {qi,q2) of the three-dimensional 
Cartesian coordinates {qi,q2,(lo), with the origin Og in the figure, under the constraint 

<il + <il + ql = \ (7) 

corresponds to two points of the sphere, where A is the curvature of the sphere. The 
points on the sphere can also be described by the spherical polar coordinate (i?, x, 0) 
defined by {qi,q2,qo) = (-Rsinxcos^, i?sinxsin6', i?cosx) with R = 1/a/A being a 
constant, (ii) The Cartesian coordinates {xi,X2) of the two-dimensional gnomonic 
projection, which is the projection onto the tangent plane from the center of the sphere 
in the embedding space, is given by 

where = xf + x\ and the point of tangency Ox in the figure being the origin. 
And the polar coordinate (r, ^) of the projection is defined by r = Rianx ^"^^ 
{X11X2) = {r cos 6, r sin 6). In this work, we mainly adopt the two coordinate systems, 
{xi, X2) and (r, 9), considering the results of Higgs [9j introduced in the following. 

In 1979, Higgs [9] introduced a generalization of the hydrogen atom and harmonic 
oscillator in a spherical space. He demonstrated that, in the gnomonic projection as 
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shown in Fig. [l] (a), the orbits of the motion on a sphere can be described by 

lL^[ry^/ + r-'] + V{r) = E-lxL^ (9) 

where the angular momentum L = X1P2—X2P1 is an invariant quantity with the potential 
V{r) being radial symmetric The Hamiltonian can be written as 

H=^^ + \\L' + V{rl (10) 

where vf = p + | [x{x ■ p) + {p- x)x] is the conserved vector in free particle motion on 
the sphere. Since the curvature appears only in the right combination E — \\L'^ of 
Eq. ([9]), the projected orbits are the same, for a given V{r)^ as in Euclidean geometry. 
Consequently, according with the Bertrand Theorem [TOlIll], the orbits are closed only 
if the potential takes the Coulomb or isotropic oscillator form, i.e. V{r) = — ^ or 
V{r) = ^u'^r'^, with k, and u being constants. Therefore the systems described by Eq. 



(10) with the two mentioned potentials are defined as the Kepler problem and isotropic 
oscillator in a spherical geometry in The algebraic relations of their conserved 
quantities reveal the dynamical symmetries of the two systems are described by the 
5*0(3) and SU{2) Lie groups respectively. These results are the beginning of the so 
called Higgs Algebra, which has been studied in a variety of directions [T2l[T3t fTH flSl fT6] . 

The concept of symmetry is one of the cornerstones in the modern physic, and 
dynamical symmetry plays a important role in many important physical models. Since 
the dynamical symmetries of the Kepler problem and isotropic oscillator on a 2-sphere 



described by Eq. (10) adhere to the behaviors in two-dimensional Euclidean geometry, 
our question is: Do there exist more qualities of being homogeneous? This paper is 
aimed at constructing the VT and the HVT for the spherical geometry and studying 
their applications. In this work, we focus on the two- and one-demential cases for 
simplicity. On the other hand, the motion on of a charged particle on a 2-sphere is not 
trivial, which is related with the famous fractionally quantized Hall states [T71IIH]- We 
provide a general equation relating the average of the kinetic energy with that of the 
potential energy in the spherical geometry. We also give a generalized HVHF, which 
could propose to solve a class of problems the sense of perturbation. 

The article is organized as follows: In the Sec. [2| the VT in both classical mechanics 
and quantum mechanics is constructed. In the Sec. |3} we generalize the VT to HVT, 
and give the quantum hypervirial relation. Two examples are taken to demonstrate the 
perturbation method which is combined HVT with HP theorem in the Sec. |4j We end 
this paper with some relevant discussions in the last section. 

2. Virial Theorem 

2.1. Classical Mechanics 



To obtain the classical VT in the two-dimensional spherical geometry, two we special 
orbits are listed in the following for examples, (i) The first one is the uniform circular 
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motion with r = as shown by the curve si in Fig. [T] (a). The kinetic energy of this 
case is given by 

T = ^{RsmxofO', 

where R sin xo is the radius of the path. The corresponding centripetal force is 

2T (l + Xr^)f-VV 



R sin Xo R sin Xo 

Hence, one can obtain 

2(r) = ((1 + Ar2)f- VV), (12) 

which can be considered as the VT under the case of uniform circular motion, (ii) The 
orbit S2 in Fig. [l] (b) depletes the case which the angular momentum L is zero. In the 
same way, the relationship between kinetic energy and potential energy can be obtain 

as 

2((1 + Ar2)r) = ((1 + Ar2)f- VV^). (13) 
These serve a good inspiration for us to presume that the VT in a spherical geometry is 

2((1 + Xr^)Tr) + 2{Te) = ((1 + \r^)r ■ WV), (14) 
where and Tq are the radial and rotational kinetic energy. 



In Appendix A, we give a proof that Eq. (14) is satisfied for an arbitrary orbit in 
the spherical space, and it equivalents to 

((1 + Ar2)f ■ VV^) = ((1 + Ar2)7r2) = ((1 + \r^){2T - \L^)). (15) 

where L = Xip2 — X2P1 is the angular momentum. It is easy to find that, when the 
curvature A — )■ 0, the above result reduces to Eq. ([T]). 

2.2. Quantum Mechanics 

In the literature [9j, to construct the the conserved quantities on the sphere, Higgs 
replaced the momentum p in the generators on the plan by the vector vr. This enlightens 
us on the subject that we can replace r-p'm Eq.(|3| by r ■ vr + vr ■ r to obtain the VT on 
the sphere. The expected value of the commutator is 

([f-7f + 7f -f,//]) =0. (16) 

For the system in the one-dimensional curve, whose Hamiltonian is given by 
H = tt'^ /2 + V with TT = p + \{x'^p + px^)/2, the above relation leads to 

((l+Ax2)^+^(l+Ax2))+^(A(l+Ax2)(l+3Ax2)) = ((1+Ax2)x^).(17) 
2 2 2 ax 



And in the two-dimensional case, from Eq.(lO) and Eq.(16), we obtain 



((l+Ar2)^ + ^(l+Ar2)) + i(A(l+Ar2)(2+3Ar2)) = {{l+\r^)r-VV).{lS) 
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In the polar coordinate, the Hamiltonian (10) can be written as 
Ho = Tr + Te + V, 



Tr 
Ta 



15 



3A+-AV2 + 
4 

1^ A — 



:i + Ar2)(l + 5Ar2) d 
r dr 



Ar 



2\2 



(19) 



where and denote the radial and rotational kinetic energy. The relation of Eq. 



(18) equivalents to 



((1 + \r^)Tr + T,(l + Ar^)) + 2(T,) + -(A(l + Ar2)(2 + 3Ar^)) = ((1 + \r^)f- VV).(20) 

These results have the same form with the classical mechanical counterparts in Eqs. 
(14) and (15), but with a term |(A(1 + Ar^)(2 + 3Ar^)) in addition. And, the term 



is different from the corresponding one in the one-dimensional case (17), which comes 
from the commutation relation of r and vf . 



3. Hypervirial Theorems 

In the above, we have got the VT in both classical mechanics and quantum mechanics. 
We will discuss the quantum HVT in the present part. A natural candidate of the 
hypervirial operator is r'^vr + vrr'^ according to r^p in the plane we mentioned in Sec. [l| 
with k being integers. 



3.1. One- dimensional 

In the one-dimensional case, one can calculate directly the commutation relation in the 
expected value 

([a;V + 7ra;^i7]) = 0, (21) 

and obtain 

[xV + 7^x^ H] = ^^{2(1 + Ax2)x^"~V2 + 2ti\1 + Ax')x'=-' 

+ (1 + \x^)[{k - l){k - 2)x'=-3 + 2AA;V-^ + \^{k + l){k + 2)x'+^]} 

-2ix\l + \x^)'^. (22) 
dx 

Because of 7r^/2 = H — V and {H) = En being the eigenvalues of the eigenstate, the 
Eq. (21 ) turns to 



2kEn{x'-')x - 2k{x'-'V)x - {x''^)x 



dx 



k, 



+ '-[{k + l){k + 2)A2(x'=+1)a + 2e\{x''-')x + {k- l){k - 2)(x'=-=^)a] = 0, (23) 

in which we denote {f)\ = ((1 + \x^)f). Hence, we get the recurrence formula of {x^)\, 
which is the quantum hypervirial relations in the one-dimensional sphere. 
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3.2. Two-dimensional 

We now consider HVT in the two-dimensional spherical geometry. For a radial potential 
V = V{r) in the Hamiltonian (19), the eigenfunction of energy can be written as 



\I^(r, 



Amd 



(24) 



with m = 0, ±1,±2... is the eigenvalue of the conserved angular momentum L. The 
Schrodinger equation 

Ho^{r,e) = E^{r,e) (25) 
reduces to the radial equation as 

Hi^{r) = Etp{r), (26) 
where the Hamiltonian Hi is given by 



d^ (l + Ar2)(l + 5Ar2) d 1 



dr 



- — + 3A + — AV^ 



It can be written as 



Hi 



+ ^1, 



+r.(27) 



f28) 



where the radial component of n is vr,. = — ^[(1 + Ar^)^ + ^ + l-^'"] ~ 
^ ~ ^[(i ~ m^)A — ~2^'^ ]- Choosing the hypervirial operator as r'^Tr,. + TTj-r^, one 
can get the recurrence relation 

2kE„{r'^~'), - 2k{r'-'Vi), - {r'^)x 



dx 



k 



+ ^[{k + l){k + 2)X^r'+')^ + 2eX{r'-')^ + {k - l){k - 2){r''-^)^] = 0, (29) 



from 



([rV + v^,r^i^l]) = 0. (30) 

Here, the notation {f)\ = ((1 + Xr'^)f). It is the two-dimensional quantum hypervirial 
relation we will discuss in the present work. And when A — 0, it reduces to the result 
in the 2-plane case [19]. 



4. Application of The Hypervirial Theorems 

In this section, we will generalize the HVHF theorem to the spherical space based on 
the hypervirial relations in the above. When the perturbation of potential V{r) takes 
the form as r'(l -|- Ar^) with / being integers, we can determine the eigenenergies in the 
various orders of approximation without calculating the wavefunction, as the the HVHF 
theorem in the Euclidean geometry. In the following, we will give two sample examples 
to illustrate this method. 
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4.I. One- dimensional Harmonic Oscillator 

The Hamiltonian of the one-dimensional harmonic oscillator in the spherical geometry 
with a perturbation potential is 

^2 1 



TC 



H =— + ^ax^ + I3x\l + Ax^), 



(31) 



where a, (3 are real numbers, / is an integer and A is the curvature of the sphere. The 
perturbation /9a;' (1 + \x^) has to be very small, and /3 is the smallness parameter. 



Then, the HVHF recurrence relation in Eq. (23) becomes 



(A; + l)a-^(fc + l)(fc + 2)A2 



k 



{x''+')x = 2kE^{x^-')^ + -X{x'^~')x + -{k - l){k - 2){x''-^)^ 

- (3{2k + /)(x'^+'-1)a - l3X{2k + / + 2)(a;''+'+i)A. (32) 

The above equation establishes precisely regarding the n-th energy level. In order to 
obtain the approximate solution of the energy eigenvalues we expand both En and 
desired expectation values {x^)\ in powers of the perturbation parameter (3 as 



En = E^^ + + + --- = Y: P'E^\ 



(33) 



3=0 



j=0 



where we introduce the notation Qj = {x^)x,j for convenience. We now insert the series 
in (33) into (32) and order in power of /3. It is straightforward to get the relation 



ik + l)a-^{k + l){k + 2)\^ 



Q 



k+l 



(34) 



—AO 
2 ^ 



fc-i 



2ky^ E^Q''-^- 

j=0 

In addition, by the HF theorem, we know that 



{k -l){k- 2)2:^-3 - {2k + l)Q';t[-' - X{2k + / + 2)Q^+[+\ 



dEn _ , dH _ I 
which gives another relationship of the coefficient of /9: 



3 



(35) 



(36) 



In other words, the j-th approximate of energy eigenvalue En^ is determined by the 
(j — l)-th approximate of desired values Qj„i. 

In the following, we would like to give an explicit example. We let / = 1 in the Eqs. 



(34) and (36) and obtain, respectively. 



{k + l)a-^{k + l){k + 2)\^ 



Q 



k+l 



j=0 



2 



2kY,EiQ'--] + ^AQ^-i + ^{k - l)ik - 2)Q^-' 



(37) 



(2A; + l)Q^_i-A(2fc + 3)Q^;?, 
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n ■ ^7 — 1 

3 

One can start from 

{x')^ = (1 + \x^) = 1 + \{x^) 

to obtain Q^. By the HF theorem 
dEn _ ,dH^, _ 1 , 2\ 



one can find that 



.dE, 



j=0 



da 



(38) 
(39) 
(40) 
(41) 



Substituting it to Eq. (39), the expectation value (x°)a expansion will be denoted as 



1 + 2A 



+ h /3 



da da da 

Ordering in power of /3, it is easy to find the first term of the recursion: 



(42) 



dE, 



(0) 



da 



1 + 



(2n + 1)A 
VA2 + 4a^ 



Q? = 2A 



Ql = 2\ 



dEiP 

da 
dE^ 
da 



(43) 



The eigenenergy of one- dimensional harmonic oscillator in a spherical geometry is 



E, 



(0) 



(n 



When 7 = 0, one can substitute Qq into Eq. (37) and obtain the values of Qq, 

k = Ql = 0, (44) 
[{2n + 1)A + VA2 + 4a] 



A; = 1 



(4a-3A2)VA2 + 4a 



(2n + 1) VA2 + 4a + (2n^ + 2n + 3)A , (45) 



Using Eq. (44) and (38), we can get the first-order perturbation of En, 
E^^ = Ql = 0. 



And from the Eq. (43) and (46), we have 

dEk'^ 



Qi = 2A- 



da 



0. 



(46) 



(47) 
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In the case of 7 = 1, using Q° Eq. (37), we can derive the values of Q{, and 



consequently the second approximation of energy level 
v/A2 + 4a + (2n + l)A 



2aVA2 + 4a 
3A[(2n + l)A + VA2 + 4a] 



(2n + l)VA2 + 4a + (2n2 + 2n + 3)A 



(4J 



2a(4a-3A2)v/A2 + 4a 
In this way, we can obtain the expectation value expansions and the energy values 
Ei^^ in the various orders of approximation as 



E^^^ = 0, 



(49) 



4a 



6A^(2^f + A) \ 

4a-3A +^^ + ^^" 



5A 



3a - 6A2 



(4E(°)+4A)Q}-(5 + 



21AM°) 
2a - 15A2 



go 

8a - 60A ° 



■(,50) 



In the limit A — > 0, En^ is tending to — l/(2a) and the other En"* is tending to zero 
which are corresponded with the exact result in the Euclidean space. 

It is worth to mentioned that, alien from Euclidean space, (i) The HP theorem has 
been used twice in this HVHF perturbative method, (ii) Only when the exponent / in 



the perturbation potential is a positive integer, we can get Qi^ from Eqs. (43) (37) and 



(38). 



4-2. Two-dimensional Coulomb System 

Here we wish to show that the HVHF perturbation method can be easily applied to 
treat the Coulomb system with a perturbation in the two-dimensional sphere which is 
described by the Hamiltonian 



TT 



K 



H=Y + 2^^' ~'r^ ^""^^ + ^^^^ 
where k is a real number, and /3 is the perturbation parameter. Hence, the potential in 



the radial Hamiltonian (28) is 



- + Ml + AO-i[(i-m^)A-^^^^] 



(52) 



The hypervirial relation Eq. (29) turns to 
1 



Xk 



^ [k(k - l){k -2)-{k- l)(4m2 - 1)] (r'^-^);, + —{k'' + 2 - Am^){r''-^)x 



+ 2kE^{r'-'), + 2{k - l)^{r'-'), + -{k + l){k + 2){r^^')^ 
- (3(2k + l){r^^^-^)x - l3Xi2k + I + 2)(r'=+'+i)A = 0. 



(53) 
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Considering the angular quantum number as a parameter of the potential Vi, one 
can obtain the expansion coefficients for {r~'^)x by using the HF theorem, 



(54) 



From this starting point, as we show in the one-dimensional case, we can get any order 
perturbation on the energy level, with the precondition that / is a negative integer. 
Taking / = —3 for example, in the first approximation, the eigenvalue En is 



Er, 



K 



1 — (n + 



^rn?){n + vm? + -) 



^m?{4m'^ — l){n + y/m? + 



/3- 



2kX 



'm2(4m2 — 1) 



(4n + 47^+1). (55) 



When A — i- 0, this result is coincided with the literature 



5. Conclusion and Discussion 



The VT in a spherical geometry has been proved in both classical and quantum 
conditions. We also have considered the HVT and got the hypervirial relations. The 
HVT and HF theorems have been shown to provide a powerful method of generating 
perturbation expansions. We have taken the Coulomb problem and harmonic oscillator 
for instances to illustrate this method. When the curvature A is zero, the results reduce 
to the counterpart of Euclidean space. 

In this paper, we only give attention to one- and two-dimensional systems. Since 
the Higgs' results have extended to the A^-dimensional spherical geometry directly [23] , 
we can foretell our treatment can be generalized to the A^-sphere and suggest the VT is 
given by {{l + Xr^)^ + ^{l + Xr^)) + l{X{l + Xr^){N + 3Xr^)) = n{{l + Xr^)f-VV). Some 
researchers have discussed the superintegrable potentials in the the hyperbolic plane [21] , 
it is interesting and possible to study the VT , HVT and HVHF in the situation of the 
curvature A < 0. On the other hand, the systems in the curved space we investigate 
in this work also can be considered as the problems with position-dependent effective 
mass, which are widely applied in various areas of material science and condensed matter 
[251 1211 [261 127] . We hope to find the applications of our results in these directions in 
the further research. 
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Appendix A. Proof of the Virial Theorem in Classical Mechanics 



In this part, we will give the strict proof for the classical VT in Eqs. (14) and (15). 
We adopt the subscripts p and s to distinguish the systems on a plane and on a sphere 
respectively. From the Eq. (|9]), we know that, for a given V{r), when 

Es — 2^L^s = Ep, Ls = Lp, (A.l) 

the projected orbit of a spherical system is the same as the orbit of a system in Euclidean 
geometry. It is easy to find that, for the corresponding points (r^, Og) = (rp, 9p) = (r, 9), 
the velocities satisfy 

Vs = il + Xr^)vp, (A.2) 

where Vs = {rs,rs9s) and Vp = {rp,rp9p). For the system in a fiat space whose 
Hamiltonian is given by H = p'^/2 + V, the two terms in Eq. ([!]) are 

{fp . vy) = - / Vp ■ VVdtp = - [fp- Vv\vp ■ dsp, 

Jo '^P Jc 

(Tp) = ^ [ \vldtp = ^ l\vl-kvp ■ dsp, (A.3) 

'p Jo ^ T^p Jc > 

where dsp = {drp,rpd9p), c denotes the orbit of motion, and Tp is the period (for the 
aperiodic case Tp — +oo). Suppose the period of the system with the same orbit c in 
the sphere described by Eq. (10) is r^. Then, considering the relations in Eqs. (A.l) 
and (A.2), one can find 

{fp-VV) = -{{l + Xrl)fs-VV), 

Tp 

{Tp) = ^[((1 + \rl)T,r) + {T,e)] = -((1 + Ar,^)^), (A.4) 

Tp Tp A 

where the radial kinetic energy Tg^ = R'^xj./2 = rf/[2(l + Ar^)^] and the rotational 
kinetic energy Tse = i?^ sin^ = r'^^l/[2{l + Ar^)]. Therefore, the relation in Eq. 

(14) is the VT in a spherical geometry, and it equivalents to Eq. (15). Here the proof 
comes to an end. 



